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We report the results of Monte Carlo simulations on several spin glass models at low temperatures. 
By using the parallel tempering (Exchange Monte Carlo) technique we are able to equilibrate down 
to low temperatures, for moderate sizes, and hence the data should not be affected by critical 
fluctuations. Our results for short range models are consistent with a picture proposed earlier that 
there are large scale excitations which cost only a finite energy in the thermodynamic limit, and 
these excitations have a surface whose fractal dimension is less than the space dimension. For the 
infinite range Viana-Bray model, our results obtained for a similar number of spins are consistent 
with standard replica symmetry breaking. 



PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 
I. INTRODUCTION 

There has recently been a renewed interest in the na- 
ture of the spin glass phase. Two principal theories 
have been investigated: the "droplet model" proposed 
by Fisher and HusetJ, (see also Refs. ||,||), and,the replica 
symmetry breaking (RSB) picture of Parisiocl. These 
scenarios and some,-pthers have also been considered by 
Newman and SteinO. An important difference between 
these models concerns the number of large-scale, low en- 
ergy excitations. RSB theory follows the exact solution of 
the infinite range Sherrington-Kirkpatrick (SK) model in 
predicting that there are excitations which involve turn- 
ing over a finite fraction of the spins and which cost only 
a finite amount of energy in the thermodynamic limit. In 
addition, the surface of these large, finite energy excita- 
tions is expectecti to be space filling which means that 
the fractal dimension of their surface, d s , is equal to the 
space dimension, d. The droplet theory argues that the 
lowest energy excitation involving a given spin and which 
has linear spatial extent L typically costs an energy L , 
where 9 is a (positive) exponent. Hence, in the thermo- 
dynamic limit, excitations which flip a finite fraction of 
the spins cost an infinite amount of energy The droplet 
theory also predicts that d s < d. n I— . 

Recently, Krzakala and Martina (KM) , and two of usE£l 
(PY), have argued that a straightforward interpretation 
of their numerical results at zero temperature is inter- 
mediate between the droplet and RSB pictures in that 
there appear to be large scale excitations whose energy 
does not increase with size, but these have a surface with 
d s < d. This interpretation of the results of KM and PY 
has, however, been recently challenged (though in oppo- 
site senses) by Marinari and ParistJj aad by MiddletonEj. 
There has also been recent debateEjO as to whether the 
± J model has similar behavior to that of a model with a 
continuous distribution. In the scenario of KM and PY 
it is necessary to introduce two exponents which describe 
the growth of the energy of an excitation of scale L: (i) 
6 (> 0) such that L is the typical change in energy when 
the boundary conditions are changed, for example from 
periodic to anti-periodic, and (ii) which characterizes 



the energy of clusters excited within the system for a fixed 
set of boundary conditions. In this paper we test whether 
the picture proposed by KM and PY is compatible with 
finite temperature Monte Carlo simulations. — __ 

Several previous Monte Carlo simulationsEZl~OEl have 
found evidence for finite energy large scale excitations 
by looking at the order parameter function P(q). In the 
thermodynamic limit this has delta functions at (plus 
or minus) the Edwards- Anderson order parameter, qEA, 
corresponding to ordering within a single valley, and, ac- 
cording to RSB theory, a tail with a finite weight extend- 
ing down to q = 0. In the droplet theory, P(q) is trivial, 
i.e. has only delta functions at ±qea, though in a finite 
system there is arweight at the origin which vanishes with 
increasing L likeEJ L~ e . 

These earlier Monte Carlo studies have found that the 
weight of P(q) at the origin is independent of the sys- 
tems size fat temperatures down to T ~ 0.7Tq-in three 
dimensionstij and T ~ 0.6 T c in four dimcjasionstJ. How- 
ever, these studies have been criticizedEj'Ej as being too 
close to the critical point, so that the results are af- 
fected by critical fluctuations and very much larger sizes 
would be needed at these temperatures to see the asymp- 
totic behavior of the low-temperature spin-glass stateE3. 
~~ HJP also argue, however, 
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that clear evidence for 
droplet theory behavior could be seen even for quite small 
sizes at very low temperatures. 

In this paper we check this prediction by performing 
Monte Carlo simulations in the low temperature region, 
though with an admittedly modest range of sizes—using 
the "parallel tempering" Monte Carlo methocoo, also 
known as Exchange Monte Carlo. One difficulty with this 
approach is to ensure equilibration since .-the technique 
proposed earlier by one of us and BhattEj for conven- 
tional Monte Carlo does not work for parallel tempering. 
Here we use an alternative method, valid for the impor- 
tant case of a Gaussian distribution (which we use here) , 
and which is closely related to the approach of Ref. |2fJ 
for the SK model. 

Both in three and four dimensions, we find a tail 
in P{q) which is independent of size (up to the sizes 
studied), for temperatures down to T ~ 0.2 T c in 3D 
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and T ~ 0.1 T c in 4D, in contrast to the prediction of 
Refs. ^0[[2lJ We also find that data for the "link over- 
lap" , defined below, fits well a description with d s < d, 
though the extrapolation to the thermodynamic limit is 
quite large here. Thus our results are completely consis- 
tent with the earlier proposal of KM and PY. 

We consider the short-range Ising spin glass in three, 
and four dimensions, and, in addition, the Viana-BrayO 
model. The latter is infinite range but with a finite aver- 
age coordination number z, and is expected to show RSB 
behavior. All these models have a finite transition tem- 
perature. In the 3D case, the exponent 8 obtained from 
the magnitude of the change of the ground state energy 
when the boundaia^conditions are changed from periodic 
to anti-periodic isc3 about 0.2, whereas in 4D it is much 
largercJ, about 0.7. 

The Hamiltonian is given by 

h= — jjjSjSj, (i) 

where, for the short range case, the sites i lie on a sim- 
ple cubic lattice in dimension d = 3 or 4 with N = L d 
sites (L < 8 in 3D, L < 5 in 4D), Si = ±1, and the 
Jij are nearest-neighbor interactions chosen according to 
a Gaussian distribution with zero mean and standard 
deviation unity. Periodic boundary conditions are ap- 
plied. For the Viana-Bray model each spin is connected 
with z = 6 other spins on average chosen randomly (but 
with the constraint that the total number of bonds is ex- 
actly 3N). We allowed the local coordination to fluctu- 
ate which is different from the more familiar Viana-Bray 
model in which each site has exactly the same coordi- 
nation number, but we expect the properties of the two 
models to be very similar. The width of the Gaussian dis- 
tribution is again unity, and the range of sizes is N < 700. 

Our attention will focus primarily on two quantities: 
the spin overlap, q, defined by 

i=i 

where "(1)" and "(2)" refer to two copies (replicas) of 
the system with identical bonds, and the link overlap, qi, 
defined by 

In the last equation, Nb is the number of bonds (Nz/2 
for the models considered here, where z is the coordina- 
tion number) , and the sum is over all pairs of spins i and 
j which are connected by bonds. The advantage of cal- 
culating qi as well as q is that if two spin configurations 
differ by flipping a large cluster then q differs from unity 
by an amount proportional to the volume of the cluster 
while qi differs from unity by an amount proportional to 
the surface of the cluster. 



II. EQUILIBRATION 

Simulations of spin glasses at low temperatures are now 
possible, at least for modest gizes , using the parallel tem- 
pering Monte Carlo methodcScJ. In this technique, one 
simulates several identical replicas of the system at dif- 
ferent temperatures, and, in addition to the usual local 
moves, one performs global moves in which the temper- 
atures of two replicas (with adjacent temperatures) are 
exchanged. It turns out to be straightforward to design 
an algorithm which satisfies the detailed balance condi- 
tion, and it will also have a good acceptance ratio if the 
temperatures are fairly close together. In this way, the 
temperature of a given replica wanders up and down in 
a random manner, and each time the temperature goes 
low the system is likely to end up in a different valley of 
the energy landscape. Thus different valleys are sampled 
in much less time than it would take for the system to 
fluctuate between valleys if the temperature stayed fixed. 

We choose a set of temperatures Ti,i — 1, 2, • • • , Nt, 
in order that the acceptance ratio for the global moves 
is satisfactory, typically greater than about 0.3. Since, 
at each temperature, we need two copies of the system 
to calculate q and qi as shown in Eqs. (||) and (||), we 
actually run 2 sets of Nt replicas and perform the global 
moves independently in each of these two sets. 

To believe results of simulations carried out at low tem- 
peratures it is essential to have a sound criterion for equi- 
libration. The technique pioneered by Ref. ^5] does not 
work with parallel tempering Monte Carlo because the 
temperature does not stay constant. However, another 
method can be used for a Gaussian distribution of ex- 
change interactions, which is a common and convenient 
choice. It depends on anJdentity first noted a long time 
ago by Bray and MooreEHl for the SK model. Here we 
give the corresponding result for the short range case. 
We start with the expression for the average energy per 
site, 

^ = -^E[ J «w^u ( 4 ) 

(*,i> 

where (• • -)t denotes the Monte Carlo average for a given 
set of bonds, and [• • -] av denotes an average over the 
(Gaussian) bonds Jij. One can perform an integration 
by parts over the J^to relate U to the average link over- 
lap defined in Eq. (||), i.e. 

where the brackets (• ■ •) indicate both a Monte Carlo av- 
erage and and an average over disorder, J 2 is the variance 
of the interactions (set equal to unity in this paper) , the 
sum is over sites i and j connected by bonds (each pair 
counted once), and the factor of z/2 arises because there 
are z/2 times as many bonds as sites. A very similar 
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approach has been used to test equilibration of the par- 
allel tempering method for the SK modeO, except that 
in that case the square of the spin glass order parameter 
appears rather than qi. 

We start the simulation by randomly choosing the 
spins in the 2Nt replicas to be uncorrelated with each 
other. This means that the two sides of Eq. (^) should 
approach the equilibrium value from opposite directions 
for the following reason. The data for (qi) will be too 
small if the system is not equilibrated because the ran- 
dom start means that the spins are initially further away 
from each other in configuration space than they will be 
in equilibrium, whereas initially the energy will not be 
as negative as in equilibrium so the right hand side of 
Eq. (g) will initially be too high. Hence we expect that 
if the two sides of Eq. (^) agree then the system is in 
equilibrium. 

For illustration purposes we show in Fig. [j] how this 
works for three-dimensions with L = 8 and T = 0.5 (to 
be compared withEl T c ~ 0.95). The data for (qi) in- 
creases as the length of the simulation increases while 
that determined from the energy decreases (by a lesser 
amount). Once the two agree they do not appear to 
change at longer times, indicating that they have reached 
equilibrium, as expected. Furthermore, the data for dif- 
ferent moments of the spin overlap q as well as the link 
overlap qi appear to saturate when qi has equilibrated so 
it does not appear that there are longer relaxation times 
for q than for qi. We also checked that the whole dis- 
tribution P(q) does not change with time once the two 
estimates for (qi) agree. By presenting results for an in- 
termediate temperature, rather than the lowest tempera- 
ture, in Fig. |l|, we show that the results do not change for 
times longer than that needed for qi and 1 — 2T\U\/z to 
agree. The length of the simulation was chosen so that, 
at lowest temperature, the data for qi and 1 — 2T\U\/z 
just converged. 

III. RESULTS 

A. Three Dimensions 

In Table [|, we show N sampi the number of samples, 
N swe ep, the total number of sweeps performed by each 
set of spins (replicas), and Nt, the number of tempera- 
ture values, used in the 3D simulations. For each size, the 
largest temperature is 2.0 and the lowest temperature is 
0.1. However, for L — 8, the data for the two lowest tem- 
peratures, T = 0.10 and 0.15, are not fully equilibrated 
so data at these temperatures has been ignored and the 
lowest temperature-used in the analysis is 0.20. This is to 
be compared withEil T c w 0.95. The set of temperatures 
is determined by requiring that the acceptance ratio for 
global moves is satisfactory for the largest size, L = 8, 
and for simplicity the same temperatures are also used 
for the smaller sizes. For L — 3, 4, 5, and 6, the accep- 
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FIG. 1. The solid squares are results for the average link 
overlap, defined by Eq. (^|), as a function of the number of 
Monte Carlo sweeps, N swccp , that each of the 2Nt replicas 
performs. Averaging was performed over the last half of the 
sweeps indicated. The triangles are obtained from the energy 
in the way indicated, and should agree with the results for 
qi if the system is in equilibrium, as shown in Eq. (^). The 
two sets of data approach each other from opposite directions 
and then do not appear to change at larger number of sweeps, 
indicating that they have equilibrated. We also show data for 
higher moments of q and qi. They appear to be independent 
of the number of sweeps once the qi data has equilibrated. 
The data for the different moments has been shifted upwards 
by the following amounts for better viewing: {q 2 } by 0.11, 
(<? 4 ) by 0.48, (qf) by 0.17, (qf) by 0.375. These results are for 
d — 3, T — 0.5 and L = 8, and the data is averaged over 3891 
samples. 

tance ratio for global moves is greater than about 0.6 in 
average and is always greater than 0.3 for each pair of 
temperatures. For L — 8 the average acceptance ratio is 
0.41, and the lowest value is 0.12. 

In Table || we compare the average total energy (E) — 
NU at T — 0.20 with the average ground state energy 
obtained by finding the ground state of each sample with 
a hybrid genetic algorithm, as discussed elsewhereEE The 
two energies are very close together, indicating that at 
this temperature our data are unlikely to be affected by 
the critical point. 

Figs. H and || show (symmetrized) data for P{q) at 
temperatures 0.20 and 0.50. There is clearly a peak for 
large q and a tail down to q = 0. At both temperatures 
one sees that the tail in the distribution is essentially 
independent of size. A more precise determination of 
the size dependence of -P(O) is shown in Fig. El where, to 
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FIG. 2. Data for the overlap distribution P(q) in 3D at 
T — 0.20. Note that the vertical scale is logarithmic to better 
make visible both the peak at large q and the tail down to 
q — 0. In this and other similar figures in the paper, we only 
display some of the data points as symbols, for clarity, but 
the lines connect all the data points. This accounts for the 
curvature in some of the lines in between neighboring symbols. 
In this paper all distributions are normalized so that the area 
shown under the curve is unity. 




samp 

15000 
16000 
7590 
4539 
3891 



■^Vj 3 WCCp 

w 

io 4 

10 5 
3 x 10 5 
10 6 



N T 
~T8 

18 
18 
18 
18 



TABLE I. Parameters of the simulations in three dimen- 
sions. iV sam p is the number of samples (i.e. sets of bonds), 
A^sweep is the total number of sweeps simulated for each of the 
2Nt replicas for a single sample, and Nt is the number of 
temperatures used in the parallel tempering method. 



L 


(E) 


(Eo) 


No 


4 


-106.49 ±0.05 


-106.60 ±0.03 


50000 


6 


-364.06 ± 0.17 


-364.94 ± 0.06 


39246 


8 


-867.04 ±0.27 


-868.20 ±0.15 


13302 



FIG. 3. Same as for Fig. H but at T = 0.50. 



TABLE II. Average energy (E) at T — 0.2 and average 
ground state energy (-Eo), for several sizes in three dimensions. 
iVo is the number of samples used to compute the average {Eo) 
using a hybrid genetic algorithm. 



improve statistics, we average over the (discrete) q- values 
with \q\ < q , with q Q = 0.20. If general, we expect that 
P(0) ~ L~ e , where we allow 8' to be different from 9, the 
latter being obtained-from boundary condition changes. 
In the droplet picturaj 8' = 8. The dashed line in Fig. || 
has slope —0.20 corresponding to the estimated valueEj of 
—8, so in the droplet picture the data is expected to follow 
a track parallel to this line. The actual size dependence 
is clearly much weaker than this, and consistent with a 
constant P(0), which implies that the energy to create a 
large excitation does not increase with size, and therefore 
8' ~ 0. 

More precisely, a two-parameter fit of the data in Fig. || 
with the form aL- ' , gives 8' = 0.01 ± 0.05 for T = 0.20, 
8' = 0.02 ± 0.02 for T = 0.34, and 8' = -0.01 ± 0.02, 
for T = 0.50. Different values of q give similar results 
for the fits. We also tried a one-parameter fit in which 
8' is fixed. Assuming 9' = the goodness-of-fit param- 
eter Q is 0.34,0.77 and 0.57 for T = 0.20,0.34 and 0.50 
respectively, whereas assuming 9' = 0.20 the goodness- 
of-fit parameters are very low, 1.4 x 10 -4 , 2.8 x 10~ 6 and 
7.6 x 10~ 15 respectively. Hence, just considering statis- 
tical errors for the sizes studied, the data is compatible 
with 9' = and not with 8' = 0.20. 

Refs. [20|j2]] studied P(0) by the Migdal-Kadanoff ap- 
proximation, which is known to yield the droplet picture 
asymptotically. They find that, although the behavior of 
P(0) at higher temperatures is masked by critical point 
effects, data at low temperatures, such as those con- 
sidered here, should show the droplet behavior. That 
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we find quite different results indicates that the Migdal- 
Kadanoff approximation is not applicable to such small 
sizes. However, our data still do not rule out the possibil- 
ity that the droplet theory, or some other theory, might 
be correct at larger sizes. 
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FIG. 5. Data for P(0) as a function of temperature in 3D 
for different values of L. The inset shows P(0)/T vs. T. 



FIG. 4. Log-log plot of P (0), the spin overlap at q — 0, 
against L in 3D. The data is independent of size within the 
error bars. The dashed line has slope —0.20, which is the 
estimated value of —6. Asymptotically, the data should be 
parallel to this line according to the droplet theory. 

In Fig. U we show the data for P(0) versus T. We see 
an approximately linear decrease of the data as T — ► 0. 
Note though, that there is some non-linearities as shown 
in the figure's inset. 

Fig. g shows the distribution of the link overlap qi at 
T = 0.20. We see that there is a large peak at qi close to 
unity (with structure coming from the allowed discrete 
values of qi) and a much weaker peak (note the loga- 
rithmic vertical scale) for smaller qi which grows slowly 
with increasing L and moves to larger values of q\. We 
will refer to this feature again below when we discuss the 
Viana-Bray model. 

The variance of P(qi) is shown in Fig. ||for several low 
temperatures. The data is consistent with a power law 
decrease to zero, i.e. 

Var(qi)~L-"«, (6) 

where the power m seems to vary somewhat with T. This 
variation is probably due to corrections to scaling coming 
from the shift with T of the complicated peak structure 
at large q seen in Fig. ||. 

The asymptotic value of \xi is related tor-the exponents 
9' and d s that have been mentioned earlierS To see this, 
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FIG. 6. The distribution of the link overlap in 3D at 
T = 0.20 for different sizes. Note the logarithmic vertical 
scale. 
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FIG. 7. Same as for Fig. g but at T = 0.50. 

assume that the non-zero variance is largely due to the 
excitation of a single large cluster of size of order L. Its 
energy is of order L e , and the probability that thermal 
fluctuations can create it is of order T/L e , assuming a 
constant density of states for these excitations. One mi- 
nus the link overlap between the two states is of order 
j^-(d-d s ) Decause there is only a contribution to 1 — qi 
from the surface of the cluster. Hence there is a prob- 
ability TL~ e of getting a Sqi of order L~^ d ~ ds \ so the 
variance goes like TL~^ 1 where 



+ 2(d~d s ). 



(7) 



An extrapolation of our results for fii to T — gives 
0.76 ± 0.03. Assuming 9' — this implies 



d- d. 



0.38 ±0.02 



(8) 



which is consistent with the T = result of PY, namely 
d-d s = 0.42 ±0.02. 

We have also looked at more general fits of the form 



FIG. 8. Log-log plot of the variance of qi as a function 
of size in 3D at several temperatures. The data for T = 0.50 
and T = 0.70 are multiplied by 1.2 and 1.7 respectively for 
better viewing. The data for T = 0.70 is somewhat curved 
and so does not fit well a power law. 



One sees that the minimum of x 2 is for a = and 
the range of a in which x 2 has increased by less than 
unity relative to the a = value is a < 5.3 x 10~ 4 for 
T = 0.50 and a < 1.3 x 10~ 3 for T = 0.34. The width of 
the distribution of qi is \fa which has values 0.023 and 
0.036 respectively. Thus while, our data cannot rule out 
a non-zero value for the width of the distribution of qi in 
the thermodynamic limit, it does suggest that this value, 
if non-zero, must be very small. We note, however, that 
a rather small value of a is not unreasonable in RSB. For 
example, if P{qi) consists of two delta functions at a dis- 
tance of 0.1, whose weights are 0.1 and 0.9 respectively, 
then the value of a is 0.0009. 



B. Four Dimensions 



Var(qj) = a + 



(9) 



to see to what extent the data can rule out a non-zero 
value of a. This is of interest because a non-zero value 
for a is required by standard replica symmetry breaking 
theory. We carried out fits of the form 



ln[Var(qi) — a] = lnfe — clnL, 



(10) 



in which a is fixed, and In b and c are the fit parameters. 
The x 2 °f the fit is then determined as a function of a 
and the results are shown in Fig. 0. 



In four dimensions we presertU,results down to a tem- 
perature of 0.20, compared withE-3 T c w 1.80. Parameters 
of the simulations are shown in Table III. For each size, 



the largest temperature is 2.80 and the lowest is 0.20. 
The acceptance ratio for global moves is always greater 
than about 0.5 for L = 3, about 0.2 for L = 4 and about 
0.3 for L = 5. In Table [Fv| we compare the average en- 
ergy at T = 0.2 with the average ground state energy 
obtained with the hybrid genetic algorithm of Ref. |3^. 
As in 3D, the data indicate that the system is very close 
to the ground state. 
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FIG. 9. The x 2 of the fit in Eq. @ (in which a is fixed 
and In 6 and c are fit parameters) for different values of a. 
The results are for the 3D model. 



FIG. 10. Data for the overlap distribution P(q) in 4D at 
T = 0.20. The data is normalized so the area under the curve 
is unity. Hence P(q) is half as big as it it would be if we had 
just plotted the region of positive q as in Figs. H and H. 



L 


N M 

1 v samp 


3 WCCp 


N T 


3 


60000 


6 x 10 3 


12 


4 


30000 


6 x 10 4 


12 


5 


12190 


3 x 10 5 


23 



TABLE III. Parameters of the simulations in four dimen- 
sions. W Quantities involving the link overlap qi have been 
calculated with half the number of samples. 



L 


(E) 


(Eo) 


No 


3 


-158.64 ±0.03 


-158.82 ± 0.04 


45000 


4 


-510.42 ± 0.07 


-510.93 ± 0.08 


26681 


5 


-1253.41 ±0.20 


-1254.33 ±0.20 


8990 



TABLE IV. Average energy (E) at T — 0.2 and average 
ground state energy (-Eo), for several sizes in four dimensions. 
No is the number of samples used to compute the average 
(Eo) using a hybrid genetic algorithm. 




FIG. 11. Same as for Fig. but at T = 0.46. 
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Figs. |h] and [fl] show data for P(q) for temperatures 
0.20 and 0.46. As in three-dimensions, the tail in the 
distribution is essentially independent of size. We display 
the full P{q), rather than just the symmetric part as in 
3D, in order to show that it has a symmetric form as 
expected (a symmetric form was also obtained in 3D). 
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FIG. 13. Data for P(0) as a function of temperature in 
4D for L = 3,4, 5. The inset shows P(0)/T vs. T. 



FIG. 12. Log-log plot of P(0) against L in 4D averaged 
over the range \q\ < 0.20. The data is independent of size 
within the error bars. The dashed line has slope —0.70, which 
is the estimated value of — 0. Asymptotically, the data should 
be parallel to this line according to the droplet theory. 

The size dependence of P(0), averaged over the range 
\q\ < q with q = 0.2 is shown in Fig. The dashed 
line has slope —0.70 corresponding to the estimated value 
of —9. In the droplet picture, the behavior should follow 
this form asymptotically. Clearly it does not for this 
small range of sizes. More precisely, performing a similar 
analysis as in 3D we find 9' = 0.10 ±0.12, 0.08 ±0.09 and 
0.17 ± 0.06 for T = 0.20, 0.32 and 0.46 respectively. For 
the same temperatures, the goodness-of-fit parameter is 
0.67, 0.65 and 0.014, assuming 9' = 0, which is accept- 
able, while assuming 9' = 0.70 the goodness-of-fit param- 
eters are tiny: 2.1 x 10- 10 ,2.1 x 10~ 12 and 4.2 x 10~ 18 . 

As in 3D we find different results from what is pre- 
dicted by Refs. ^0|,^l] at such low temperatures on the 
basis of the the Migdal-Kadanoff approximation. How- 
ever, our data can not rule out the possibility that some 
other behavior may occur at larger sizes. 

Note also, that our data for P(0) decreases approxi- 
mately linearly with temperature as T — > 0, as shown in 
Figure 

Figs. |14| and |l5| show the distribution of the link over- 
lap qi at temperatures 0.20 and 0.46. As in 3D we see 
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FIG. 14. The distribution of the link overlap in 4D at 
T = 0.20 for different sizes. Note the logarithmic vertical 
scale. 
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FIG. 15. Same as for Fig. |lj but at temperature 0.46. 

complicated structure at large qi and a subsidiary peak 
at smaller qi which grows with increasing L. 
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temperatures. The data is consistent with the power law 
decrease to zero shown in Eq. (||). The range of sizes 
is so small, and the values of \xi also so small, that we 
are not able to rule out a non-zero value for L — ► oo in 
4D. However, the data is consistent with the asymptotic 
value being zero. 

An extrapolation of our effective values of fii to T = 
gives 0.35 ± 0.06, which, assuming 0' = 0, gives 

d-d s =0.17 ±0.03. (11) 

This is just consistent with the T = results of PY 
who find d — d s = 0.21 ± 0.01. However, the quoted error 
bars are from statistical errors only, so the difference may 
be partly due to systematic effects coming from the small 
range of sizes studied. 

C. Viana-Bray Model 

For the Viana-Bray model, T c is given by the solution 

of 

-L f°° e~ x2 ' 2 tanh 2 ( -) dx = — , (12) 

where z (= 6 here) is the coordination number. The 
solution is T c = 1.8075 which is roughly twice the 
transition temperature of the 3D short range model con- 
sidered here, which has the same coordination number. 

Parameters of the simulations are shown in Table [v|. 
In each case, the largest temperature is 2.6 and the lowest 
temperature is 0.1. For N = 700 the data is not equili- 
brated for temperatures lower than 0.34, and is almost 
equilibrated at T = 0.34. Except for N = 700 the accep- 
tance ratio for global moves is always greater than about 
0.3. For N — 700 the acceptance ratio is greater than 
0.3 for most temperatures but there is one "bottleneck" 
where the acceptance ratio went down to 0.08. 



N -/Vsamp ^Vswcep 

59 19022 IF 21 

99 5326 3 x 10 4 21 

199 3116 10 5 21 

399 3320 10 5 21 

700 801 3 x 10 5 21 



TABLE V. Parameters of the simulations for the 
Viana-Bray model. 



FIG. 16. Log-log plot of the variance of qi as a function 
of size in 4D at several temperatures. The data for T = 0.63 
and T — 0.81 are multiplied by 1.15 and 1.7 respectively for 
better viewing. 

The variance of qi is shown in Fig. ftu at several low 



First of all, in Fig. [i~7| we show that P(q) has a weight 
at q = which appears to be independent of the system 
size, as expected. 

A plot of P(qi) is shown in Figs. [TJ and |l| for T = 0.34 
and 0.70. Note that the data for N = 700, T = 0.34 show 
a dip around q ~ 0.5, due to imperfect equilibration of 
the Monte Carlo runs for some samples (this explains 







also the fluctuations of the data for P(q) in Fig. 17). 
As in the 3D and 4D data discussed above we see a 2- 
peak structure develop as the size increases. This is also 
clearly visible in the earlier work of Ciria et al.E!l in 4D. 
For the Viana-Bray model, the position of the smaller 
peak shifts neither with L or T whereas for the 3D data, 
see Fig. [j| it clearly shifts to larger values of qi with 
increasing size. In 4D, see Figs. and the range of 
sizes is sufficiently small that it is difficult to tell whether 
there is a shift in the position of the peak or not, but if 
it is present, it appears to be a smaller effect than in 3D. 
It would be helpful to understand the physics behind the 
two peak structure. For the hierarchical lattice r»sed in 
the Migdal-Kadanoff approximation, Bokil et al£H have 
given an explanation, but it is not clear to us how this 
goes over to the models discussed here. 




FIG. 17. Data for the distribution of the overlap for the 
Viana-Bray model at T = 0.34. Note the logarithmic vertical 
scale. 

A plot of Var(qi) against N is shown in Fig. ^o] for 
different temperatures. In contrast to the data for 3D 
shown in Fig. |^, (which is for a similar number of spins) 
the data is clearly tending to a constant at large N 

This is confirmed by the x 
responding to Eq. ( |l0| ) shown 
asymptotic value of a is large and finite. Compare this 
figure with Fig. |9| which shows the corresponding results 
in 3D, and where a = gives the best fit. 

We conclude this section by pointing out the sizes stud- 
ied (which covers a similar range to that in 3D and 4D) 
are sufficient to determine the correct asymptotic behav- 
ior for the Viana-Bray model. 



2 analysis of the fits cor- 
in Fig. |l|. Clearly the 




FIG. 18. Data for the distribution of the link overlap for 
the Viana-Bray model at T = 0.34. Note the logarithmic 
vertical scale. 



Viana-Bray, T = 0.70 




FIG. 19. Same as for Fig. h8 but for T = 0.70. 
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IV. CONCLUSIONS 
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FIG. 20. A log-log plot of the variance of qi for the 
Viana-Bray model at different temperatures. The data at 
T = 0.5 and T — 0.7 are multiplied by 1.2 and 1.6 respec- 
tively for better viewing. 
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FIG. 21. The x 2 of the fit in Eq. (0) for different values 
of a for the Viana-Bray model. Note that the minima of the 
different x 2 are not monotonic in T since we have fewer data 
points for T = 0.20. 



To conclude, Monte Carlo simulations at low (but fi- 
nite) temperatures agree with earlier T = studies of 
KM and PY that there appear to be large-scale low en- 
ergy excitations which cost a finite energy, and whose 
surface has fractal dimension less than d. However, since 
the sizes that we study are quite small, there could be 
a crossover at larger sizes to different behavior, such as 
the droplet theory (with 9' = 9 (> 0)) or an RSB picture 
(where 9' = 0, d — d s = 0). We note, however, that our 
results for short range models are quite different from 
those of the mean-field like Viana-Bray model for sam- 
ples with a similar number of spins, and, furthermore, our 
results for the Viana-Bray model do predict the correct 
asymptotic behavior for that model. 
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